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Hydrostatic equilibrium of a porous intracluster medium: 
implications for mass fraction and X-ray luminosity 
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ABSTRACT 

The presence of dilute hot cavities in the intracluster medium (ICM) at the cores 
of clusters of galaxies changes the relation between gas temperature and its X-ray 
emission properties. Using the hydrostatic equations of a porous medium we solve for 
the ICM density for a given temperature as a function of the filling factor of dilute 
bubbles. We find that at a given temperature, the core X-ray luminosity increases with 
the filling factor. If the frequency of AGNs in clusters were higher in the past, then 
the filling factor could correspondingly be significant, with implications for the cluster 
scaling relations at high redshifts. This is especially important for the core properties, 
including the temperature-luminosity [Lx — T) relation and estimates of the core gas 
mass. The results imply an epoch-dependent sensitivity of the Lx — T relation in the 
core to the porosity of the ICM. Detection of such an effect would give new insights 
into AGN feedback. 
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1 INTRODUCTION 

X-ray observations of the intracluster medium (ICM) of 
massive clusters indicate the presence of cavities seen as de- 
pressions in the X-ray surface brightness (e.g. Boehringer 
et. al. 1993; Fabian ei. al. 2000; McNamara et. al. 2000; 
Birzan et. al. 2004). These cavities or bubbles are believed 
to have been generated by AGN activity of central galaxies 
and possibly other galaxies in the cluster (Nusser, Silk & 
Babul 2006; Best et al. 2007; Martini, Mulchaey, & Kelson 
2007). 

The observed properties of the x-ray emitting gas in 
clusters of galaxies provide strong constraints on the thermo- 
dynamic properties of the ICM. The most satisfactory model 
which accounts for the observed self-similar scaling proper- 
ties of outer regions invokes galaxy formation by cooling gas 
(Suginohara & Ostriker 1998; Lewis et al. 2000; Pearce et al. 
2000; Dave, Katz & Weinberg 2002; Ettori et al. 2004; Nagai, 
Kravtsov & Vikhlinin 2007), but this in turn leads to a major 
puzzle. Approximately twice as much cold gas is predicted as 
compared with the mass in stars. AGN feedback in the core 
can resolve this problem (e.g. Quilis et. al. 2001, Babul et. 
al. 2002, Dalla Vecchia et. al. 2004, Ruszkowski et. al. 2004, 
Volt & Donahue 2005, Briiggen et. al. 2005, Roychowdhury 



et. al. 2005, Sijacki & Springel 2006, Nusser, Silk & Babul 
2006, Peterson & Fabian 2006, McNamara & Nulsen 2007, 
Cattaneo & Teyssier 2007). Related discussions of the role 
of AGN heating of the ICM are given in semi-analytic mod- 
elling of galaxy formation (Croton et al. 2006; Bower et al. 
2006). We point out here that AGN feedback drives poros- 
ity in the ICM and modifies the hydrostatic properties of 
cluster cores. A higher AGN activity in the past leads to a 
redshift evolution of the scaling relations. 



2 THE HYDROSTATIC EQUATION FOR THE 
AMBIENT MEDIUM 

A formalism for following the evolution of bubbles statisti- 
cally has been developed in Nusser, Silk & Babul (2006). 
In this formalism, the hot dilute bubbles and the rest of the 
ICM are a two-fluid system that is best described as a single 
fluid which we term the ambient medium. The representa- 
tion in terms of a single medium is applicable only if local 
pressure equilibrium between the bubbles and the residual 
ICM is established. Hereafter, quantities related to the ambi- 
ent medium, the bubbles, and the residual ICM (gas outside 
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bubbles) are denoted by the subscripts, a, b, and /, respec- 
tively. Suppose that the volume filling factor of bubbles at 
r is F(r). The ambient density, p^, at position r, is defined 



Pa = (1 - F)p, + Fp, 



(1) 



in terms of the density inside bubbles, Pj^ , and the density 
of the ICM, Pj, also at r. The local ambient energy, u^, per 
unit mass as 



(1 - F)p,u, + Fp, 

Pa 



(2) 



where Uj^ and Uj are the energy per unit mass of the material 
inside the bubbles and of residual the ICM. If the pressure, 
p, is p = (7 — l)itp for both the ICM and the material inside 
the bubbles. The equation of state of the ambient medium 

is 



P = (7a " l)Pa"a > 



where 



7a 



(7i - 1)(7. - 1) 



(1-F)(7, -1)+F(7, -1) 



(3) 



(4) 







(5) 



and we have assumed local pressure equilibrium between 
the bubbles and the ICM. Hydrostatic equilibrium of the 
ambient medium in a gravitational field g{r) is described 
by, 

1 dp 
Pa dr 

We numerically integrate these equations assuming 
NFW profile with a virial radius of 4 Mpc and a concen- 
tration parameter of 7 (Neto et al. 2007). We assume that 
the gas outside the bubbles has a temperature of 6 x lO'^K 
and that the core baryonic to dark mass ratio equals 0.14. 
We take a core radius, re = 100 kpc. This choice of the 
parameters yield reasonable temperature and density in the 
core. Of course, the choice is not unique. Once can obtain 
similar characteristic of the central region for smaller virial 
radii at the expense of changing the concentration and the 
core radius. For simplicity we assume that the mass in the 
bubbles is negligible. 



3 RESULTS 

The density averaged over the bubbles and residual core 
gas in the ambient medium, p^, as a function of radius is 
plotted in figure (1) for various values of the filling factor. In 
each case, a constant filling factor is assumed. The profiles 
become flatter as the fllling factor is increased. This can 
be seen analytically from the hydrostatic equation (5) as 
follows. Substituting p^ = (1 — F)p-^ and P = (7 — l)PiWi in 
this equation we obtain, for it, = const and F = const, 



P = (7 - 1)W] 



dlnpj 



1/(1-F) 



dr 



(6) 



This implies that the density Pj [F) as a function of the filling 
factor scales like Pi{F) --^ Pi{F — 0)^~^ which in turn leads 
to a similar scaling for the ambient density plotted in the 



figure. Since < F < 1, this scaling yields flatter profiles 
for larger F. 

The modification to the density profile causes a 
change in the corresponding core X-ray luminosity, Lx oc 
/ drr^(l - F)Pi^(r) = / drr^(l - Fy^pJ{r). The X-ray 
luminosity (normalized to the value for F = 0) is plotted in 
figure (2). The luminosity is an increasing function of the 
filling factor. This is due essentially to the larger ambient 
density in the outer regions for a porous medium, amplified 
by the (1 — F)^^ factor in the expression for Lx- The net 
result is a boost in the core X-ray luminosity. 



4 DISCUSSION 

Our results may be interpreted as follows: 

(i) For a given observed X-ray emission and spectroscopic 
temperature of the residual ICM gas (outside the bubbles) , 
hydrostatic equilibrium can yield erroneous mass estimates 
if the porosity of the medium is not properly included. 

(ii) The Lx ~ T relation will depend on the porosity. The 
effect is enhanced if analysis is restricted to the core where 
the porosity is highest. This is in addition to changes re- 
sulting from gas cooling in the core (Allen & Fabian 1998, 
Markevitch 1998). 

(iii) The ICM porosity is expected to be larger in the past 
due to an enhanced frequency of AGN activity. This leads 
to redshift dependence of a core Lx — T relation which is 
different from that seen for the outer parts of the cluster 
(Kotov & Vikhlinin 2005). 

(iv) Porosity enhances the X-ray emission and hence 
shortens the cooling time. This leads to more AGN activity 
and the effect is only quenched when the porosity becomes 
sufficiently large, when only little amount of gas remains 
available for feeding the AGN. At this point the cold gas is 
all driven out of the core until there is time for a new episode 
of cooling and associated AGN activity. 

Predictions of the modified hydrostatic equation should 
be compared against observations by the analysis of hydro- 
static equilibrium of clusters. Porosity of individual clusters 
is needed for such a study. But reliable estimates of the 
porosity is not always possible using current data. One could 
alternatively proceed by solving the modified equation for 
assumed porosity as a function of cluster temperature and 
mass. The scatter in the relevant relations could then be 
minimised with respect to the assumed porosity values. 

Rough estimates of porosity could be derived from the 
results of Birzan et al. (2004) (see also Dunn & Fabian 2004 
and Dunn, Fabian & Taylor 2005, Dunn & Fabian 2006). 
Using their table 2 for the sizes of prominent active bubbles 
(radio lobes), we compute the filling factors corresponding 
to some of their clusters. The fllling factors within a distance 
of 30 and 40kpc from the center are listed in table 1 here. 
The range of the filling factors corresponds to the uncer- 
tainty in projection effects (cf. Bfzan et al.). Each cluster in 
table 1 contains 2 radio lobes and the upper (lower) limiting 
value of the fllling factor is obtained by summing the max- 
imum (minimum) values of the sizes of both lobes. Table 1 



© 0000 RAS, MNRAS 000, 000-000 



A porous intracluster medium 3 













j;^o 


4.6 


- 


^^~~~--~-„„^^^ 






F=0.2 ^~~~~~-~.^^ 






F=0^6^^^^-C^\\ 


4.4 


- 


F=0.8 ^~~~~~~-~-~.,^\V\ 






~~~~~"^^S5^ 


4.2 


_ 


^ 









-0.8 



-0.6 -0.4 

log(r/rJ 



-0.2 



Figure 1. The ambient density, p^, profiles (in units of the baclcground value) as a function of radius (in units of re, for several values 
of the bubble filling factor. 




Figure 2. The normalized X-ray luminosity versus the filling factor. 



© 0000 RAS, MNRAS 000, 000-000 



Nusser & Silk 



Table 1. Filling factors of radio lobes in a few rich clusters (data taken from Birzan et al. 2004). 



F(_R < 30kpc) 
¥{R < 40kpc) 



Perseus 
0.006 -0.036 
0.002-0.015 



Hydra A 
0.07-0.14 
0.06-0.11 



Cygnus A 

0.11-0.25 

0.12-0.26 



MKW 3S 
0.12-0.23 
0.09-0.19 



RBS 797 
0.023 -0.14 
0.01 -0.066 



A2199 

0.006-0.03 

0.002-0.01 



demonstrates that the porosity resulting from active promi- 
nent bubbles could be significant. Note that, low contrast 
bubbles remain undetected (Dunn & Fabian 2006) so that 
the porosity could be much higher than the estimate based 
on observed bubbles. Therefore, more rigorous estimates of 
the porosity, especially at larger radii and at different red- 
shifts, by detailed observational studies of higher S/N are 
needed in order to test the evolutionary implications of the 
AGN feedback hypothesis. 

The modified hydrostatic equation is valid for any two- 
component fiuid where the density of one of the components 
is negligible. The generalization to multi-component fiuids 
of any densities could easily be done using the derivation 
described in Nusser, Silk & Babul (2006). 
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